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Abstract. We have developed a Monte Carlo simulation describing the 2^2 
scatterings of perturbatively produced, non-eikonally propagating high-energy partons 
with the quarks and gluons of the expanding QCD medium created in ultrarelativistic 
heavy ion collisions. The partonic scattering rates are computed in leading-order 
perturbative QCD (pQCD), while three different hydrodynamical scenarios are used to 
model the strongly interacting medium. We compare our results and tune the model 
with the neutral pion suppression observed in ^snn = 200 GeV Au+Au collisions at 
the BNL-RHIC. We find the incoherent nature of elastic energy loss incompatible with 
the measured angular dependence of the suppression. The effects of the initial state 
density fluctuations of the bulk medium are found to be small. Also the extrapolation 
from RHIC to the LHC is discussed. 



1. Introduction 

Our Monte Carlo (MC) simulation [TJ [2J [3] aims to study the pQCD-based collisional 
energy loss of hard partons traversing a strongly interacting medium. It is similar 
to the MC models JEWEL (Jet Evolution With Energy Loss) [4] and MARTINI 
(Modular Algorithm for Relativistic Treatment of heavy IoN Interactions) [5], but while 
JEWEL and MARTINI include both elastic and radiative energy-loss components, we 
concentrate purely on the energy loss induced by the elastic 2 — >■ 2 scattering processes. 
In this study, we aim to answer the following questions: i) Within reasonable parameter 
values, how large could the elastic energy loss contribution be? ii) Does the pathlength 
dependence of elastic (or, in general, incoherent) energy loss match the experimental 
data? hi) What are the effects of initial state density fluctuations to energy loss? 

2. The model 

We model the elastic energy loss of a hard parton by incoherent partonic 2 — > 2 processes 
in pQCD with scattering partners from the medium, upgrading our hydro as we move 
to the more complex observables. For central heavy ion collisions, a (l+l)-dimensional 
hydro [6] with initial conditions from the EKRT model [7] is used to describe the 
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medium. For non-central collisions, we utilize a (2+1) -dimensional hydro [H] with a 
smooth sWN profile j9] obtained from the optical Glauber model. Finally, to study 
the effects of the initial state density fluctuations, we have an event-by-event hydro [8] 
with an eBC profile [9] from the Monte Carlo Glauber model. Additionally, we use 
(2+l)-dimensional hydro with a smooth eBC profile pU] for LHC simulations. 

The basis of our approach is the scattering rate r i (p 1 ,u{x),T(x)) for a high- 
energy parton of a type i with 4- momentum pi, accounting for all possible partonic 
processes ij — > kl. The hydrodynamical model provides the local flow 4-velocity u(x) 
and temperature T(x) of the medium. In the local rest-frame of the fluid, we can express 
the scattering rate for a process ij — > kl as follows PQ: 

Here E\ is the energy of the high-energy parton i in this frame and E2 is the energy of 
the thermal particle j with a distribution function fj{E2, T), which is assumed to be the 
Bose-Einstein distribution for gluons and the Fermi-Dirac distribution for quarks. The 
scattering amplitude \M\f^ kl depends on the standard Mandelstam variables s,t,u. A 
thermal-mass-like overall cut-off scale m = s m g s T is introduced in order to regularize 
the t, u -singularities appearing in the scattering amplitude. The mass factor s m and 
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the strong coupling constant a s = j*- are the free parameters of our model. We keep a s 
fixed with momentum scale. 

The hard parton is propagated through the plasma in small time steps At. At 
each step, the probability for a collision is given by the Poisson distribution 1 — e~ TiAt . 
After scattering, the parton with higher energy is selected as the hard parton. As we 
assume no significant interactions between the high-energy parton and the fully hadronic 
medium, no collisions are allowed in regions with temperature below the decoupling 
temperature Td ec - The medium-modified distribution of high-energy partons obtained 
in the end can be convoluted with a fragmentation function to calculate the standard 
nuclear modification factor R AA (P T , y , <p) = {N ^^) d Zvv/dP Tdyd4> for a given centrality. 



3. Results 



We focus first on the high-F T neutral pions produced in ^/snn = 200 GeV Au+Au 
collisions. The simulations for 0-10% central collisions with the (l + l)-dimensional 
hydro background show that roughly the right amount of nuclear modification in the 
most central collisions can be achieved by setting a s = 0.5 and s m = 1.0 (see also the 
left panel of Figure 1). We keep these values as we proceed to non-central collisions 
with the (2+l)-dimensional hydro to examine the pathlength dependence. 

As seen from Figure [1] however, our model cannot reproduce the reaction-plane 
angle dependence seen in the PHENIX experiment. Also, the computed suppression 
decreases too slowly as one advances to the more peripheral collisions, which means 
that the inclusive, angle-averaged nuclear modification factor does not match with the 
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Figure 1. (Color online. Originally published in [2]-) Left panel: The tt° nuclear 
modification factor for 0-10% centrality, averaged over the reaction plane angle. Middle 
and right panel: The 7r° nuclear modification factor dependence on the reaction 
plane angle A(f> for 40-50% (middle panel) and 50-60% centrality (right panel). The 
simulation points (solid squares and triangles) are connected with lines to guide the 
eye. The PHENIX data are from [11] (0-10% centrality, open squares) and [12] (40- 
50% and 50-60% centrality, open triangles). Colored bars with small cap represent 
statistical errors; black bars with wide cap are systematic errors. 



experimental data either. Thus the pathlength-dependent observables strongly disfavor 
any large contributions from these kind of fully incoherent processes to the energy loss 
of hard partons. Utilizing our event-by-event hydro, we have also studied the Raa for 
fluctuating initial state geometry, where the weak sensitivity to the angle-dependent 
observables is also seen. Overall, our findings, combined with the similar results for 
radiative energy loss, suggest that the initial state density fluctuations do not play a 
major role in hard parton energy loss [3]. 

Figure [2] shows how the Raa(Pt) changes as we move from ^/snn = 200 GeV 
Au+Au collisions at RHIC to ^/snn = 2.76 TeV Pb+Pb collisions at the LHC. Using a 
(2+l)-dimensional hydro with a smooth eBC profile (as opposed to sWN profile) requires 
us to re-tune a s to 0.55 in our model to get the right amount of suppression for neutral 
pion Raa(Pt) in RHIC energies. Given the large systematic uncertainties of the present 
data, the Pt dependence of our elastic energy loss model is in reasonable agreement 
with both the RHIC and published LHC measurements. However, comparison with 
the preliminary data [13] suggests that the obtained Pt slope at the LHC may be too 
modest. Fixing the pathlength dependence without changing the Pt dependence too 
dramatically will be an interesting challenge for future studies. 
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Figure 2. (Color online) Comparison of simulation and experimental data for neutral 
pion Raa(Pt) m RHIC ^/snn = 200 GeV Au+Au collisions and charged hadron 
Raa(Pt) in LHC ^/snn = 2.76 TeV Pb+Pb collisions. The simulation points (solid 
squares and triangles) are connected with lines to guide the eye. The PHENIX data 
(open triangles) are from pj] and ALICE data (open squares) are from [14] . Error bars 
with small cap represent statistical errors; bars with wide cap are systematic errors. 
The black histogram indicates the lower alternative pp reference for the ALICE data. 
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